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We have carried out theoretical studies of the phonons and spinons excitation spectra of 
Sr2CuOs in term of a one-dimensional model with non-adiabatic interaction magnetic and 
elastic subsystems using continuous time Monte Carlo method. Two gaps in triplet excitation 
spectrum are found. The bond defects concentration estimated is well agreement with the 
experimental data. The formation of magnetic moment on site at low temperature Tjy ~ 
5K is caused by freezing quantum spin liquid with correlation radius of antiferromagnetic 
ordering £/a ~ 300 — 500. 

PACS: 75.10.Jm, 75.30.Fv, 67.20. +k, 63.20.Ry 

1. Introduction 

Quasi-one-dimensional spin systems with antiferromagnetic interactions have received con- 
siderable attention because of a number of basic physical phenomena in low dimensions, 
such as the existence of a spin gap depending on the spin value or the frustration anti- 
ferromagnetic exchange, spin Peierls transitions, accompanied by the structural symmetry 
breaking. Orthorhombic compound Sr%CuO^ is considered as quasi-lD antiferromagnetic 
system. However, Raman studies reported that extra peaks were observed along the Cu-0 
chain directions [|IJ which implies that there might be some structural instabilities. This ef- 
fect is confirmed by polarization-dependent infrared-active phonon spectra which have shown 
splitting of the stretching mode in the Cu-0 chain direction about 20cm -1 and became evi- 
dent as T decreases [|2| . These unusual behavior are not expected by the previous theoretical 
calculations. 

Another unclear behavior is the uniform spin susceptibility = 0) have decreased steeply 
below T ~ 33f^ without the signature of three -dimensional short range order approaching 
the Neel state by NMR Knight shift at the in-chain oxygen. The nuclear spin-lattice relax- 
ation rate offsets the ~ 1/T increase of the staggered spin susceptibility down to 1§K and 
depends on frequency. The q = mode of spin transport is ballistic and a diffusion-like con- 
tribution at T << J H . The bulk susceptibility the measurements have shown decreasing 
above ~ 10i^ Q . The NMR lineshape have been changed in a very specific way, giving 
rise to "features", which have been explained by the chain edges and the presence of mobile 
bond-defects. A "defect" is defined as a local change in the magnetic bond coupling. As 
the translational symmetry of the spin system is broken, a local alternating magnetization 
develops around the defects in the direction of the applied field. The density of bond-defects 
decreases linearly with T. The source of the bond-defects might be related to the interstitial 
excess oxygen. However , at the lowest temperature, T = 20K, the density of bond-defects 
p ~ 2.8 • 10~ 3 remains more than ten times larger than the concentration of the residual 
spin-1/2 impurities (~ 1.3 ■ 10~ 4 ) after annealing and the observed T dependence of p is not 
a constant . 
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The concept of mobile "bond-defects" have been used for explanation of the temperature 
dependence and the anisotropy of the thermal conductivity below 50 K as the main source 
of spinon scattering at low temperatures ||. The calculated mean distance between two 
neighboring bond defects, consistent with the NMR data between 20 and 60 K. The results of 
the specific heat measurements have been fit to the data above 4 K using the approximation 
C = 7T + (3T 3 + 5T 5 . The linear in T contribution is due to spinons excitations of a 
ID Heisenberg spin- 1/2 antiferromagnetic system with exchange J = 2620 ± 100 which is 
18% larger value than J ~ 2200 deduced from magnetic susceptibility data [|] and NMR 
data The Neel temperature Tjv dependents on the measurement time and shifted to 
lower temperature with increasing time. So Neel temperature is T/v = 5AK from neutron 
scattering |7j , 4.15K by fiSR and 3.5K from the anomalies of the specific heat 0. 
Possibly , here it is nonequilibrium process. Physical origin of the bond-defects remains 
unclear also as abrupt change suseptibility at different temperatures determined by NMR 
and by external magnetic field. 

In the paper we try to answer on these questions basing on non-adiabatic interaction between 
spins . Numerical methods, such as exact diagonalization and density-matrix renormalization 
group, face the potential difficulty in dealing with a very large Hilbert space of the phonons. 
Here , we use Monte Carlo approaches restricted to finite chains L = 100, 200 but without 
any adiabatic approximation and the truncation of the infinite phonon Hilbert space. The 
method |J is based on a path-integral representation for discrete system in which we work 
directly in the Euclidean time continuous. All the configuration update procedures contain 
no small parameters . Being based on local updates only, it allows to work with the grand 
canonical ensemble and to calculate any dynamical correlation function, expectation values. 
2. Model and method 

The ratio intechain to intrachain exchange in Sr 2 CuO^ is an extremely small ~ 10~ 3 and 
demonstrates a good one-dimensionality. The structural fluctuation were only observed along 
in chain at low temperature. As a result a spin-phonon interaction were considered along 
chain. The model Hamiltonian of an spin-phonon system takes the form: 

H = £[J + a{ Ul - u i+1 )) [StS z l+l + (S+Sr +1 + SrS+ + 1 )/2] + Mu 2 /2 + K( Ui - u l+1 ) 2 /2. (1) 

Here S z,± are a spin operator components associated with the site i, J > is the usual 
antiferromagnetic exchange integral, a is the spin-phonon coupling constant, U{ is the dis- 
placement in the z- direction, M is the mass of the ion and K the spring constant. Using 
the quantum representation for phonon operators b, b + , the Hamiltonian becomes 

# = EE[^+«y^y^ cos (^^ (2) 

+ '£ i m(q)b+b q , fi(g) = 2^sin(|), ^ = 2^, q = 2nn/L,n = 1,2, ...L 

The phonon frequency (ouq), spin-phonon coupling (a), energy E, temperature T are nor- 
malized on the exchange J and h — 1,M — 1. The temperature used in calculation is 
(3 = J/T = 25. Our system consists of the two subsystem interacted. The elastic subsystem 
is described by phonons with the number occupation n ph = 0, 1, 2.. and magnetic subsystem 
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is characterized by the number occupation n m = 0,1 and Pauli operators a, a + which coin- 
cide with S spin operators. We start with the standard Green function of the phonon in the 
momentum q- imarginary-time r representation: 

G(q, r) =< vac\b q {r)b+{Q) \vac >, r > 0, b q (r) = exp (HT)b q exp (—Ht) . (3) 
G (Q^ r ) = I < v\b}\vac > | 2 exp \-{E v {q) - E )t] 

V 

where \v > is a complete set of the Hamiltonian H in the sector of given q , H\u(q) >= 
E v (q)\v{q) >, H\vac >= E \vac >,E = 0. Rewriting Eq.(3) as 

/•CO 

G(q,r) = duA(q,uj)exp-(ujT),A(q,uj) = V<f(w - E v {q))\ < v\b+\vac > | 2 (4) 
Jo v 

one defines the spectral function A(q, uj) . We calculate a one-particle spin Green function 

< a(r)a + (0) >= / dujp t (uj) exp — (cor), (5) 
J o 

where p*(a;) -spectral density function. Correlation functions are determined on the basis of 
a complete set of eigenvector of the Hamiltonian: 



< O >= ^^'Y 1 ^ , (6) 





v 3 > 




Vj > 



where O is the longitudinal spin-spin correlation function < S^S^ +h > and the phonon 
density- density < n ph (q)n ph (q + p) >. The distribution of phonons number n p h(q) and the 
magnons number n m (k) as a function of momentum are simulated. 

3. Results and discussion 

The upper boundary of acoustic phonon band o>o is in direct proportion to Debye temperature 
of Qsr 2 CuO-i = 441±10iT in term of Debye approximation. Using known data for CuGeO^ 



, ®CuGe0 3 = 330if 13, cu 0iGe ~ 1U0K |nj we estimate u ~ 2100iT for Sr 2 Cu0 3 . The 
ratio of ujq/J ~ 0.95 will be used for calculation of dynamical characteristics . The spectral 
density of triplet excitations are shown in Fig.l. At the critical value of spin-phonon coupling 
a c i the spectral density p t = for E < A, were A is the gap in the triplet excitation 
spectrum. The gaps energies presented in inset of Fig.l fit satisfactory on the straight line 
A = —0.17(2) + 1.16(5)a in the range of a c i < a < a c2 with a cl = 0.15(1). 

For small spin-phonon coupling a < a c \ the composite quasiparticle consisted spinon and 
phonon are formed. The inelastic interaction between spinon with momentum k and phonon 
with q lead to new quasiparticle with energy E = E q +Ek and momentum k+q. It follow from 
the distribution of the phonon number versus momentum as plotted in Fig. 2a. The density of 
spinon excitations exhibits the singularity at the E(k = vr/2) and phonon excitations at the 
E(q = 7r) and summered momentum agrees with MC result N p h(q — 4.5) ^ 0. Quasiparticles 
with q — k are absent . It may be attributed to the repulsion spin-phonon quasiparticle with 
momentum q — k from spinon with the same energy E = E^ — E q . Whereas the density 
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of spinon excitations tends to zero for energies E = E q + E^. The number of spin-phonon 
quasiparticles is increased against a as illustrated in Fig. 2 and the bound states of spin- 
phonon quasiparticls are formed at a > a c \ that are similar to bipolaron with the symbolic 

kind < b+ „ ni . 6+ ™ , nr i+n2 si >. 

The correlation function of the phonon numbers < N p h(k)N p h(k + q) >, presented in Fig. 3a 
, is nonzero at the certain wavenumber q ~ 0.1 and q ~ 7r/3 for a ~ a c \ and for all q 
at the condition a > a c2 . The elastic strains cause nonuniform distribution of exchange 
and spin density that confirms spin-spin correlation function. The difference between spin- 
spin correlation function simulated in chain with spin-phonon coupling and without it is 
shown in Fig. 3b. Small oscillations are observed for a < a c2 . Cycle of the oscillations 
5r = 9 well agrees with cycle of oscillations of the magnons number distribution function 
5k = 0.36, 0.86; 5r = it /5k ~ 9 plotted in Fig. 4. The distribution of phonon numbers 
becomes continuum at the a > a c2 with one sharp maximum at the wave vector of structure 
being in the range of 7r < Q < 2ir as plotted in Fig. 2b. The value a c2 = 0.4 ± 0.02 is similar 
to a critical concentration when the dressed phonons effectively "percolate " along the three 
low-lying bands. The correlation function < N ph (k)N ph (k + q) > is not equal to zero for all 
momentums as shown in Fig. 3. 

The transition from one state with localized phonons to delocalized is like to Anderson's tran- 
sition in disordered systems. The spin density in chain allocates random and nonuniform 
that follow from the normalized spin-spin correlation function shown in Fig. 3b. The esti- 
mated values fit by exponential approximation |< SqS* > — < SqS* >af|~ Aexp(— r/£), 
where £ is a correlation radius of inhomogeneity. The magnons number distribution function 
exhibits also the random oscillations (Fig. 4c). The transition to delocalized phonons is ac- 
companied by strong increasing of the mean-square amplitude of ions vibration as presented 
in Fig. 5b. At a = a c \ the function < u 2 > (a) reveals also a some peculiarity. The average 
number of phonons , plotted in Fig. 5a , increases vs. a according to the power law 

N av = 0,0012(— ) 1 - 75 ( 6 )(7) 

«c2 



From the spectral density of phonons shown in Fig. 6 and spin excitations it may be stated 
that the coherent nonlinear bound excitations of phonon and magnon exist at a > a c \. It may 
be interpreted as a standing vibration with the spin kink in sites. The disjoint vibrations are 
performed at the condition that the wave length is changed in two times l/a = 2; 4; 8; 16; ... 
c k = n/l,E = w sin(k/2). These estimates of energy are in good agreement with Monte 
Carlo results up to l/a = 32 because the temperature fluctuations (T = 0.04) cut the kinks 
interaction radius. 

Now we will apply our results to estimate spin-phonon coupling for the one-dimensional spin 
system Sr 2 Cu0 3 with intrachange J = 2200i^, u = 2100/T. The abrupt change of the bulk 
susceptibility at T ~ 10K is associated with gap in the triplet excitation spectrum. Using 
typical ratio between gap and critical temperature, for example, for CuGeOs, A/T c ~ 1.8 
we estimate a = 0.154 from eq. 0.008 = —0.17 + 1.16a. The average phonons num- 
ber is N p h = 2.2 • 10~ 4 . These phonons are formed bound states with energy about 
~ 0.55 eV ,0.76 eV, 1 eV as followed from the spectral phonon density plotted in Fig. 6 
because maximum phonon energy with n q — 1, E — 0.95J ~ 0.18 eV. This energy of elastic 
fluctuation is compared to the value of the experimental charge transfer gap E g = 1.5±0.3eV 
derived from He I ultraviolet photoemission combined with bremsstrahlung isochromat spec- 
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troscopy [12] and the on-site Coulomb interaction on oxygen site U p = 4AeV |1| . The elastic 
fluctuation lead to change of the electron density on oxygen site and discontinuity of anti- 
ferromagnetic exchange that causes formation of the paramagnetic spin 1.3 • 10 -4 observed 
at low temperature ~ 4K. 

There are spin-phonon quasiparticles with energy E SjP h = 1040cm -1 , 562cm -1 , 287cm -1 , 

144cm -1 , 72cm -1 , 36cm -1 , 18cm -1 . The optic phonon mode observed at uj\ = 550cm -1 
corresponds to stretching mode and mode with uo 2 = 569cm -1 is attributed to spin-phonon 
quasiparticle with k = tt/4,uj mc = 562cm -1 . The relaxation time (or lifetime ) of the 
bound spin-phonon quasiparticle is decreased at heating and at T ~ E ~ 800if it is 
decomposed that qualitatively agrees with experimental data. The phonon width of optical 
conductivity spectrum is increased from 7 = 17.4cm -1 , T = 10K to 7 = 35.1cm -1 , T = 
300f^ and the strength of the oscillator is decreased from S = 0.09cm -1 eV,T = 10f^ to 
S = 0.06cm -1 el / , T = 300K. The phonon width and the strength of the oscillator of mode 
u)\ = 550cm -1 don't change in the range 10-fT < T < 300K" 0. 

The bond defects observed by NMR and heat transport arise from formation of the standing 
vibration sites. The density of bond defects p ~ 1/7 is in direct proportional to wavevector 
of the standing vibration k = 71/ 1 and estimated density is satisfactory fit with experimental 
data at T = 20K,p MC ~ 5 • 10 -3 ,p ex ~ 2.8 • 10 -3 . If discrete spectrum of spin-phonon 
quasiparticles map to continuum the temperature dependence p may be presented as p ~ 
2T / ttujq which well agree with NMR data []5| . It is also explained the linear dependence of the 
specific heat at low temperatures. Interaction between loosely coupled phonon and spinon 
modes in the center of Brillouin band may be presented as (k— u/v ph )(k — uj/v sp ) = a 2 , where 
v p h = ujo/2, v sp = n/2J, a = a/ ' ■s/uJq. It lead to repulse of these branches and formation of gap 
in the spinon spectrum with A(k = 0)/J = ay/n j2\fj ~ 0.019. This gap A(k = 0) = 43K 
causes abrupt decrease of the dynamical spin susceptibility for k = at T ~ 33f^. The 
k = mode of spin transport has a diffusion like contribution for T << J as a result of the 
fact that the finite density of bound spin-phonon quasiparticles exists at E — . 

At low temperature the spin liquid with non-linear topological excitations having S = 1/2 
is freezed that lead to small magnetic moment on site and elastic neutrons scattering. So 
in compound CuGeO^ having transition singlet-triplet at T = 14K a long range order is 
observed at low temperatures when ions Cu or Ge are displaced by Mg or Sr. Breaking 
a exchange between the nearest neighbors results in a finite spin density near breaking. 
Between the generated magnetic moments the interaction is existed per phonon field in 3D 
dimensional space. The critical field may be estimated by mean field (MF) approximation 
T N = 2Sc(J + 2Ji), where c— impurities concentration, S = 1/2 , J and Ji is accordingly 
intra- and interchain exchange in simple cubic lattice. The calculated critical temperatures 
for Cu\_ x Mg x GeO^ are in well agreement with experimental data Tff F = 2AK, T^f = 
2.5K,x = 0.016; T^ IF = 3.2K,T^ X = 3K, x = 0.0216 @ . Using similar evaluation for 



Sr2CuOs with concentration of bond defect at T ~ 6K, c ~ 2 • 10 -3 we determined the 
critical temperature Tjf ~ 4.5i^ which is in well agreement with experiment 3.5K < Tff < 
5AK. Scatter of temperatures may be due to dependence of freezing bond defects from 
frequency. The average distance between bond defects is ~ 2 • 10 3 A and the periodic local 
alternating magnetization exists which is observed by NMR at low temperatures. Mean 
square displacement is < U 2 >~ 10 -6 that may give rise the structural deformation along 
chain 5a ~ V< U 2 > ~ 10 -3 A for a = 3.9A, M ~ 10 -2 V 
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In summary, the spectrum of spinon excitations along chain is asymmetric relative wavenum- 
ber k = 7r/2 with two gaps A(k = 0) = 43K and A(k — tt) — 18K that account for steeply 
decrease of the bulk susceptibility and the dynamical spin susceptibility for k = at the 
temperatures T ~ 10K, T ~ 33K. The spectrum of bound spin-phonon excitations is due to 
bond defects the density of which is linearly increased against temperature. At low temper- 
ature T ~ 5K the non-linear spin-phonon excitations are freezed on the breaking exchange 
and induce magnetic moment on site with correlation radius of antiferromagnetic ordering 
1200 — 2000A The magnetic properties of Sr 2 CuOz are similar to CuGeOs with impu- 
rities having two transitions on temperature. The paramagnetic contribution in the bulk 
susceptibility at T < 4K result from a structural defects on oxygen ions due to spin-phonon 
coupling. 

The author are grateful to V.A. Kashurnikov, A.S. Mishchenko, N.B. Prokof'ev, B.S. Svis- 
tunov for useful assistance and discussions of Monte Carlo method in the continuous time. 
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Captions to Aplesnin paper " Quantum spin liquid in Sr 2 Cu0 3 induced by Acoustic Phonons 

» 

Fig.l The spectral density of the one-particle excitations with o> = 0.95, a = 0.05(1), 0.15(2) (a); 
a = 0.2(1), 0.4(2)(6). In insert, the gap energy A in the excitation spectrum as a function 
of spin-phonon coupling . 

Fig. 2 Distribution of phonon numbers versus momentum for a = 0.05(1), 0.2(2) (a); a = 
0.4(1), 0.6(2)(6). 

Fig. 3 The normalized correlation function of the phonon numbers versus wavenumber for 
a = 0.2(1), 0.4(2), 0.6(3)(a) and the spin-spin correlation function as a function of distance 
for a = 0.2(1), 0.6(2) (b) . Solid line is shown the exponential dependence Aexp(— r/£), A = 
0.008(8), £ = 9.1(9). 

Fig. 4 Distribution of magnon numbers versus momentum for a = 0.05,0 < k < n (a); 
a = 0.2, < k < 7T (h); a = 0.6, < k < ir (c). 

Fig. 5 The average phonons number < N p ^ av > (a) and the mean-square displacement of 
ions < u 2 > (6) versus a/a c2 . 

Fig. 6 The normalized spectral density of phonon excitations p p h/p P h,max versus energy for 
the spin-phonon coupling a = 0.25(a); a = 0.4(6); a = 0.7(c). 



8 




E 



0,0012 -, 



0,0008 - 

A 
a; 

SZ 
CL 

^ 0,0004 - 



0,0000 



* |1 I) j 







0,08 - 



0,06 - 



A 



V 



0,02 - 



0,00 




J q 

2 






b 


■ 1 




^ J 




1 1 1 1 1 1 1 1 1 




1 1 1 1 1 







q 




0,009- 



A 

N 

CO 

c 

N 

CO 

V 



A 0,006- 

N 

CO 

o 

N 

CO 

y. 0,003 -I 



o 



0,000 







■ 1 

o 2 




o 



10 



15 



20 



r 



0,25 -| 
0,20- 
0,15 
0,10- 
V 0,05- 
0,00 



2" 



o!o 0^5 iV 



— i 1 1 1 1 1 1 1 1 

1,5 2,0 2,5 3,0 3,5 



0,012- 



0,008- 



V 0,004- 



0,000 




0,0 
0,012 n 



0,008- 



A 
J*. 



0,004- 



0,000 



— i < 1 1 1 1 1 < 1 1 1 < 1 

0,5 1,0 1,5 2,0 2,5 3,0 3,5 




0,0 0^5 IX) 1^5 2^0 2^5 3^0 3,5 



0,0020 



A 



> 



V 



0,0015- 



0,0010- 



0,0005- 



0,0000 



-i — i — i — i — i — i — i — i — i — i — i — i — i 
0,2 0,4 0,6 0,8 1,0 1,2 1,4 



2,0x10 i 
1,5x10" 5 - 

A -5 

^-3 1,0x10 H 
V 

5,0x1 0" 6 H 



0,0 



a/a 



c2 



0,0 0,5 



1,0 1,5 



2,0 



a I a 



c2 



